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ABSTRACT 

We establish a simple model for the distribution of cold gas around L* galaxies using a large set of observa- 
tional constraints on the properties of strong Mg II absorber systems. Our analysis suggests that the halos of 
L* galaxies are filled with cool gaseous clouds having sizes of order 1 kpc and densities of ^ 10^^ cm^^. We 
then investigate the physical effects of cloud irradiation by a quasar and study the resulting spectral signatures. 
We show that quasar activity gives rise to (i) extended narrow-line emission on 100 kpc scales and (ii) an 
anisotropy in the properties of the absorbing gas arising from the geometry of the quasar radiation field. Pro- 
vided that quasars reside in halos several times more massive than those of L* galaxies, our model predictions 
appear to be in agreement with observations of narrow emission-line nebulae around quasars and the recent 
detections of ^ 100 kpc cold gaseous envelopes around those objects, suggesting a common origin for these 
phenomena. We discuss the implications of our results for understanding absorption systems, probing quasar 
environments at high redshifts, and testing the quasar unification scheme. 

Subject headings: atomic processes — galaxies: halos — quasars: absorption lines — quasars: emission lines 
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1. INTRODUCTION 

Observational data on quasar absorption lines indicates that 
L* galaxies are surrounded by cool 10^ K) gas condensa- 
tions ranging up to radii of ^ lOOh^Q kpc^ (Bergeron et al. 
1987, Steidel et al. 1995, Churchill et al. 2005, Zibetti et 
al. 2006). Similar distributions have recently been found 
around quasars (Bowen et al. 2006, 2007). In addition, sev- 
eral teams have discovered that bright quasars are surrounded 
by giant line-emitting nebulae observed at transitions such as 
Lya (e.g., Wampler et al. 1975, Heckman et al. 1991a,b, 
Christensen et al. 2006). In this paper we investigate whether 
the two phenomena could originate from the same type of gas 
clouds around galaxies thereby connecting gas absorption and 
line emission in galaxy and quasar halos. 

Quasar spectra have historically been a main source of in- 
formation about large scale matter distribution in the universe. 
Such information has been obtained by studying absorption 
lines induced by intervening material between us and the 
quasar. In the optical wavelength range, absorption lines by 
metal-enriched gas at z < 2 originate from low to moderate- 
ionization species and the easiest feature to detect is usually 
the Mg II doublet (AA2796. 35, 2803.53). It is the first reso- 
nance transition of an abundant element with a large oscilla- 
tor strength to enter the optical range and has therefore been 
the most commonly used tracer of cool gas. Several stud- 
ies have shown that Mg II arises in gas spanning more than 
five orders of magnitude in neut ral hydrogen column density , 
from Nm lO'^ - 10^^ cm~^ ([Be rgeron & Stasin"skal ll986l: 
ISteidel & SargenlllI992l:IChurchill e t al. 2000). Following the 
suggestion by Bahcall & Spitzer (1969) that quasar absorp- 
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tion lines might arise in halos around galaxies, an association 
between absorbers and galaxies has been reported, starting 
from small sarnples of a few to ~ 50 galaxies CBergeron 1986: 
ICristianillI987l: iBergeron & Boisselll99Tl ISteidel et alj|1997l: 
Churchill et al. 2005) to statistical measurements involving 
thousands of objects (Zibetti et al. 2005 & 2006, Menai'd et 
al. 2007). 

While a connection between galaxies and metal absorbers 
has been firmly established, the origin of the absorbing gas 
remains unclear. Whether these absorption line systems trace 
gas being accreted by a galaxy or outflowing from it is still a 
matter of debate. Among the possible origins for absorption 
line systems in galaxy halos are condensations of cool clouds 
from hot and dilute halo gas via thermal instability (e.g.. Mo 
& Miralda-Escude 1996, Mailer & Bullock 2004), cool gas 
that is bound to dark matter sub-halos embedded within the 
main halo (e.g., Sternberg, McKee, & Wolfire 2002 and refer- 
ences therein), and starburst driven winds from galaxies (e.g., 
Oppenheimer & Dave 2006, Prochter et al. 2006). Addi- 
tional explanations include (warped) galaxy disks (Bowen et 
al. 1995, Prochaska & Wolfe 1998) as well as low surface 
brightness companion objects (York et al. 1986, Petitjean & 
Bergeron 1994). 

Some galaxies contain active nuclei and it is thought that 
all quasars are harbored in galaxies. The properties of quasar 
hosts are not well characterized but it is believed that bright 
quasars are usually harbored in rather massive, star-forming 
galaxies ranging from ~ 1 to ^ 30L* (Jahnke & Wisotzki 
2003). Little is known about the environments of quasars, 
with recent studies suggesting that brightest objects live in 
somewhat over-dense environments (Serber et al. 2006). Like 
galaxies, the immediate environment of quasars can also be 
probed in absorption: using a survey of projected quasar pairs 
from the Sloan Digital Sky Survey (SDSS) our team has re- 
cently showed that quasars, like galaxies, are surrounded by 
Mgll gas on scales ranging up to IOO/z^q' kpc (Bowen et al. 
2006). Using a larger sample of quasar pairs, we have now 
found that the transverse gas distribution around quasars, as 
traced by Mg II, suggests a covering factor of order unity up 
to - 100 /i^o kpc (Bowen et al. 2007). 
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Gas which is associated with the quasar environment is also 
commonly observed (via resonance line absorption such as 
C IV A 1548) in our line-of-sight toward quasars and is identi- 
fied by its proximity (within ^ lO-' km s~^) to the quasar red- 
shift. Nevertheless, it is generally hard to establish whether 
such gas lies in the vicinity of the black hole (on sub-pc 
scales) or whether it originates in its host galaxy or its sur- 
rounding halo (see Crenshaw et al. 2003 for a review). 

Deep imaging and spectroscopic observations of quasars 
reveal that a significant fraction of them have surrounding 
emission-line nebulosities extending out to ^ IOO/z^q' kpc. 
Such nebulosities are often detected in Lya emission (e.g., 
Heckman et al. 1991, Christensen et al. 2006) or the 
[O III] A5007 Une (Stockton & MacKenty 1987, Fu & Stock- 
ton 2006) with recent surveys suggesting that radio-loud 
quasars (RLQ) may have brighter Lya nebulosities (Chris- 
tensen et al. 2006). Little is known about the physical nature 
of these nebulae and various scenarios pertaining to both gas 
accretion onto the quasar as well as gas ejection from it have 
been put forward as possible explanations (e.g., Heckman et 
al. 1991, Haiman & Rees 2001). 

Motivated by these recent findings, we investigate in this 
work whether the emitting gas seen around quasars is con- 
sistent with being the same material as that observed around 
galaxies (and quasars) in absorption. The paper is organized 
as follows: In §2 we first summarize some observational con- 
straints on the properties of cool gas around galaxies used to 
calibrate our model. A geometrical model is then constructed 
and its physics is discussed. We then (§3) use this model to 
study the effects on the absorption and emission properties 
of halo gas when exposed to a quasar radiation field. In ad- 
dition, we show that such a model can naturally explain the 
distribution of cool gas around quasars and the existence of 
large-scale Lya emitting nebulae. We discuss the implica- 
tions of our results for galaxy formation, baryon fraction in 
the Universe, and quasar physics in §4. Surmnary follows in 
§5. 

2. THE DISTRIBUTION OF COOL GAS 
AROUND L* GALAXIES 

2.1. Observational constraints from Mg II absorbers 

Observational data on Mg U absorption lines have provided 
a number of constraints on the spatial distribution and charac- 
teristics of the cool gas around galaxies. Here we list the most 
relevant results that can be used to model the Mg II distribu- 
tion in the halo of an L* galaxy: 

1. ) Strong Mg n systems, defined with a rest equivalent 
width Wo > O.3A, inhabit the halos of ^ L* galaxies, as was 
found by Steideletal. (1997). These authors found that Mg II 
selected galaxy luminosities lie in the range 0.3 — 5L* with a 
mean (L) ~ 0.8L*. Their results also suggested a unit cover- 
ing factor up to ~50 hjf^ kpc and the absence of strong Mg II 
absorption on larger scales. More recent studies (Zibetti et al. 
2005 & 2006, Churchill et al. 2005) indicate that strong Mg II 
absorption can actually be detected on scales reaching more 
than IOO/z^q' kpc but that the covering factor falls off in this 
range. We note that the study done by Zibetti et al. (2006) 
implies little evolution in the properties of absorbing galaxies 
up to z ~ 1. 

2. ) Strong Mg II absorbers usually show saturated absorp- 
tion lines (as implied by the doublet ratio). High-resolution 
spectroscopy reveals that strong systems usually break into 
several kinematically distinct components (e.g., Churchill & 



Vogt 2001). When resolved, the widths of these components 
are consistent with gas at a temperature of ^ 10^ K (e.g., 
Lanzetta & Bowen 1990). It has been estabUshed that the rest 
equivalent width of strong systems is roughly proportional to 
the number of distinct kinematical components. Bergeron & 
Petitjean (1990) found the relation between total rest equiva- 
lent width and number of clouds to be Wo (Mg II) ~ 0.3 x A'c A 
and ChurchiIl(1997)reportedWo(MgII) -0.07 xA^cAwhere 
A'c is the number of clouds along the line-of-sight. If the in- 
dividual clouds observed in these systems are relatively simi- 
lar, then these results suggest that, even if the rest equivalent 
width of Mg II systems is mostly driven by a velocity disper- 
sion, it can be used as a proxy for the number of clouds along 
the Une of sight. 

3. ) Constraints on the sizes of Mg II systems have been 
obtained by studying the spectra of strongly lensed quasars. 
Limits on the size of individual absorption components are 
obtained by comparing absorption features across different 
lines of sight. Typical values span the range 0.1 — 10 kpc for 
Mg II clouds (e.g., Rauch et al. 2002, ElHson et al. 2004) and 
perhaps somewhat larger sizes for systems detected by high 
ionization lines such as C IV A A 1 548 , 1 550. At present this is 
the only direct means for estimating the size of cool gaseous 
clouds in galaxy halos (for an indirect method relying on pho- 
toionization modeling see e.g.. Ding et al. 2003). 

4. ) Finally, HI measurements show that strong Mg II sys- 
tems with Wo ~ O.5A give rise to HI colunm densities in the 
range lO'^"^*^ cm~-^ while stronger systems with Wq 2A 
haveA'ffl ~ 10'^~^i cm'^ (Rao etal. 2006). Emission-based 
results from 21 cm surveys at z = (Zwaan et al. 2005) in- 
dicate that the contribution of galaxy disks to the observed 
HI column density distribution can be substantial within —20 
kpc for certain orientations. (See also Bowen et al. 1995 and 
Prochaska & Wolfe 1997). However, on larger scales, the con- 
tribution of halo gas clouds is in general expected to dominate 
the HI column density distribution. The lack of correlation 
between Mg II absorption properties and galaxy size or incli- 
nation further supports this assertion (Steidel et al. 2002; see 
also Kacprzak et al. 2007). 

2.2. Modeling the spatial distribution 

We now introduce a simple phenomenological model aimed 
at reproducing the main aspects of the observational con- 
straints presented above, namely the spatial extent of the gas, 
the cross-section for strong Mg II absorption and the distri- 
bution of hydrogen column densities observed around ~ L* 
galaxies. We emphasize that, in this study, we do not attempt 
to investigate the physical origin of the gas, i.e. whether it is 
bound to satellite dark matter halos, infalling or being ejected 
out of the galaxy via winds. We consider a case in which the 
total mass density (baryonic and dark matter) for such a sys- 
tem is described by an isothermal profile: 
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where a is the velocity dispersion, r^ the virial radius, and p 
is the density. The mass within a given radius r is then given 
by 

M{< r) = (2) 

and by definition, the density enclosed within the virial radius 
is 200 times the critical density. According to Fukugita & 
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Fig . 1 . — The number-count distribution of Mg II clouds along different lines-of-sight through a halo of an L* galaxy (see text). Note the increase in the number 
of intercepted cloud with decreasing impact parameter Note also that clouds are distributed up to the virial radius though the probability for intercepting them 
on large scales is small. Right: The covering factor for strong Mg II absorption around galaxies. The different solid lines correspond to model predictions of the 
covering factor for different Wq limits (denoted to the right of each curve) around an L* galaxy. As shown, the covering factor for strong, Wq > 0.5A absorption is 
below 50% for impact parameters greater than 50 kpc. Conversely, weak absorption with < 0.3A is common for impact parameters greater than 70 kpc. The 
dashed Unes are model predictions for galaxy halo models with cool gas mass 0.25 and 4 times that of an L* galaxy and for Wo > 0.5. 



Peebles (2006), the one-dimensional mass-weighted velocity 
dispersion for late type galaxies is oiatc = 140km s^' which 
corresponds to a virial mass ~ 2.5 x lO'^M© within = 
282 kpc. 

For simpUcity we assume that the cool gas around L* galax- 
ies is of the form 

Pcooi('') °= r~^- (3) 
Such a density profile arises naturally in a wide range of 
models for the origin of the absorbing gas. In particular, it 
characterizes gas which follows the dark matter distribution 
in isothermal halos. It also a property of mass conserving 
outflowing winds or infalling gas (in spherical geometry) that 
holds over a range of radii where the velocity is roughly con- 
stant. As we shall later see, deviations from this power-law 
do not significantly affect our conclusions. 

The multiple components observed in high-resolution spec- 
tra indicate that the gas does not follow an overall smooth dis- 
tribution but is distributed in discrete systems. To reproduce 
this property, we consider the cool gas to be in the form of 
clouds that are distributed within the halo following a Pois- 
son distribution with a mean number density, Wcioud following 
pcooi{r)- Considering a line-of-sight through the gaseous halo 
of a galaxy, the number of clouds intercepted along a physical 
path dl is 

^'A'cioud = icioud (r) Jr/?^ioud di, (4) 

where nR^ is the cloud cross-section for inducing Mg II ab- 
sorption, r is a vector from the halo center, and dl is measured 
along the line-of-sight. By integrating over the density profile 
of the gas we find the number of clouds along a given line-of- 
sight at some impact parameter, b, to be 

Noiondib) = ^^{ry/b) X ncioud(fe) 7iR^ b (5) 

where "^{r./b) - (l/n) f^'''' d^^-' / s/^^. 

As mentioned above, the overall rest equivalent width of 
strong Mg II systems is roughly proportional to the number of 
subcomponents. We can therefore use this empirical relation 
and approximate the rest equivalent width of a Mg II absorber 
by 

Wo^^cioudWc (6) 



with Wc = O.2A. Here we have used the mean of the propor- 
tionality factors reported by Bergeron & Petitjean (1990) and 
Churchill et al (2003) 

The final constraint needed for our model is the normal- 
ization of the overall amplitude. Steidel (1993) has shown 
that Mg II absorption with Wq > 0.5A is found within ^50 
kpc of L* galaxies with a unit covering factor. This imposes 
a constraint on our model: (Wb(f" = 50kpc)) ~ O.SA. Given 
the previous relation, it implies that, on average, about two 
Mg II clouds are intercepted by a line-of-sight with an im- 
pact parameter of 50 kpc. More recent results by Zibetti et 
al. (2005) & Churchill et al. (2005) have shown that there 
in fact exists a continuous distribution of impact parameters 
up to more than 100 kpc and the data suggest that the cov- 
ering factor is less than unity on those scales. We note that 
once the function Wo(r) has been normalized at 50 kpc, the 
cloud-distribution of Mg II systems introduced above auto- 
matically reproduces the drop of the covering factor for ab- 
sorption and allows absorbing clouds to be found on scales 
> 100 kpc, as indicated by recent studies. We note that the 
scatter in Wo(fe) °= A'cioud(^) is largely due to the relatively 
small number of clouds. Specifically, the standard deviation 
is given by [1 -|- ^/Wcioudl^) +0.75]Wc (Gehrels 1986) and is 
therefore considerable (AWo/Wo < 1). Interestingly, a large 
scatter is also seen in the data yet its magnitude is poorly con- 
strained due to small number statistics per galaxy luminosity 
bin (Steidel et al. 1993). 

We now combine aU the above properties and display the 
resulting cloud distribution in Fig. 1 . The left panel shows 
the distribution of intercepted clouds of cool gas as a function 
of impact parameter as well as the expected scatter. The right 
panel presents, for different absorption line sensitivities, the 
covering factor for absorption around L* galaxies. We note 
that the observational constraints we are using are valid on 
scales greater than ^20 kpc. This is denoted in Fig. 1 by the 
white dashed circle. On scales smaller than 20 kpc, the contri- 
bution from the disk of late-type galaxies becomes important. 
In this paper we are only interested in gas distribution in the 
halo and we will focus only on these larger scales. 
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Fig. 2. — The various components of the ionizing radiation field consid- 
ered in this work. The quasai' spectral energy distribution (solid line), the 
galaxy ionizing radiation field (dash-dotted line), and the UV/X-ray back- 
ground radiation field (dash line). The relative normalizations of the different 
components were scaled to the values that would be measured at a distance 
of 50kpc from the center of the galaxy (see text). 
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Fig. 3. — A contour plot of the Mg IIA2796 optical depth (blue con- 
tours) and the H I column density (red contours), both in logarithmic units, 
as a function of the cloud particle number density and total hydrogen column 
density. Only the UV/X-ray ionizing background radiation is assumed. Lines 
con'esponding to several cloud sizes are also drawn (see text). 



2.3. Properties ofMg II clouds 

We now investigate the possible range of physical condi- 
tions which allow the existence of Mg II clouds around a 
galaxy. Given the low gas temperatures associated with in- 
dividual Mg II absorption components, we neglect the effect 
of coUisional heating by shocks, conduction, and cosmic ray 
heating (e.g., Bergeron & Stasinka 1986, Charlton et al. 2000) 
and only consider the effects of photoionization. For clouds 
in a galactic halo, most of the ionizing flux originates from 
the meta-galactic field (Haardt & Madau 1996) and the cen- 
tral galaxy. Contributions to the mean ionizing photon flux 
from satellite galaxies and other diffuse sources can be ne- 
glected. The background UV/X-ray ionizing radiation is due 
to quasars and stars. Characterizing this spectrum requires 
knowledge of the star formation rate and quasar activation 
history as well as the opacity of intervening material. In this 
work we use the Haardt & Madau (1996) calculations and 
consider the background radiation at z '--^ 1 . The spectrum, in- 
tegrated over 4n steradians, is shown in Fig. |2]with a dashed 
line. 

In addition to diffuse background radiation there is also 
the ionizing radiation from the central galaxy. We model the 
galaxy spectrum using the Bruzual & Chariot (1996) models 
and the Chabrier (2003) initial mass function. We assume 
grey opacity beyond the Lyman edge so that 5% of the ion- 
izing continuum escapes to the halo (e.g., Bland-Hawthorn 
& Maloney 2001). A population of X-ray binaries is also 
included (Norman et al. 2004) using the infra-red to X-ray 
scaling of Ranalli et al. (2003). The galaxy ionizing radiation 
is shown in Fig. |2] 

To model the photoionization and thermal equilibrium of 
the cool gas clouds in the halo, we use the CLOUDY c06.02a 
photoionization code (Ferland et al. 1998) with the Badnell 
di-electronic recombination coefficients (Badnell 2006). We 
simplify the calculations by assuming a slab geometry with 
one surface of the cloud exposed to half of the total ionizing 
flux. Photoionization calculations are carried out to the center 
of the cloud and symmetry is assumed when calculating the 
total column density of ions (see Gnat & Sternberg 2004 for a 



more realistic treatment of the problem). We assume the gas 
composition to be 0.1 x solar (Turnshek et al. 2005, York et 
al. 2006). 

We investigated several possible scenarios for the structure 
of the clouds: (a) constant density; this is perhaps the simplest 
model yet has been the main route by which a wide range of 
astrophysical phenomena are modeled, (b) constant radiation 
and gas pressure; this may be relevant to clouds in pressure 
equilibrium with their surroundings (e.g., with the inter-cloud 
medium), (c) constant temperature; this could be envisioned 
as gas which is collisionally heated to the virial temperature of 
dark matter minihaloes. Given the parameter space of interest, 
it appears that the different models do not give appreciably 
different results. For clarity we will show only the results 
obtained for a constant density medium. 

We have calculated the column density of Mg II and H I 
as a function of cloud hydrogen (ionized and neutral) column 
density and the number density of particles. The results are 
shown in Fig. [3] 

We have repeated the calculations including the effect of 
the ionizing radiation field of the galaxy and found it to be 
dominant within the central ^ 40kpc (see §2.4 for the effect 
of cloud shielding). The gas considered here, having a low 
ionization level, is sensitive to the ionizing flux just above 
the Lyman edge. To first order, gas which is exposed to the 
galaxy ionizing flux (in addition to the UV ionizing radiation 
flux, /uv), fg would attain a similar ionization structure to 
that shown in Fig. [3]provided its density is greater by a fac- 
tor 1 +/g//uv- Similarly, at higher redshifts, where the UV 
background radiation is stronger the density of a cloud with 
a similar ionization level to that shown in Fig. |3] would be 
higher by the UV flux ratio of the epochs. 

Constraining the parameter range relevant to strong Mg II 
absorbers is not straightforward since strong Mg II lines are 
usually saturated and additional observational constraints are 
needed to estimate their typical density and column density. A 
useful constraint is the size of the clouds derived from strong 
lensing studies. As noted above, such studies indicate that 
the size of Mg II absorbers is ^Ikpc. This constraint sub- 
stantially reduces the relevant parameter space in Fig. |3] A 
second constraint comes from the H I column density mea- 
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Fig. 4. — HI column density vs. the rest equivalent width of Mg 11 for 
Mg II selected systems (see Rao et al. 2006). Data appear as black points 
with upper limits marked by triangles. Different point size indicates different 
selection criteria employed by Rao et al. (lai'ge points are systems selected 
by criteria 1-2 while smaller points by 3-4; see Rao et al. 2006); both sets of 
points have similar properties for the low column density range considered 
here. The contiibution of galaxy disks to the measured H I column density 
may become important in the shaded region of the diagram (Zwaan et al. 
2005) and our model does not attempt to account for it. The simulated data 
are drawn in red points and reproduce the lower column density envelope (see 
text). 

surements for strong absorbers (Rao et al. 2006). Using the 
rest equivalent width as a proxy for the column density along 
the line-of-sight (see eq. ^ we see from Fig. |4]that a typical 
cloud has an H I column density of lO'^ cm^^, implying a 
density of about 2 x 10^^ cm^^. Such values are consistent 
with previous estimates in the literature for the density of cool 
gas in galaxy halos (e.g., Sternberg, McKee, & Wolfire 2002 
for the case of high velocity clouds in the Galaxy). 

The scatter in the Nhi vs. Wq relation presented in Fig. |4] 
is significantly larger than that of the simplified model pre- 
sented in the previous section. Several effects not included 
in our model could also contribute to this scatter: a distribu- 
tion of galaxy masses, types, star formation histories, etc. - 
all of which are beyond the scope of this paper and are not 
observationally motivated at this stage. In the absence of a 
clear physical motivation, we attempt to take this effect into 
account by introducing a distribution of cloud sizes (as im- 
plied by the observations; §2.1). Here, clouds have similar 
densities which allow them to be in pressure equilibrium with 
their environment. (We find the alternative scenario in which 
all clouds have the same size but different densities less plau- 
sible since the implied evaporation times would be embarrass- 
ingly short compared to the Hubble time; see appendix.) To 
this end we assume that the size of individual Mg II clouds is 
drawn from an intrinsic distribution ^ (/^cioud) so that observed 
one (taking into account cloud cross-section effects on the de- 
tectability rates) 0obs(/^cioud) °^ ^(/^cioud) x ^cioud- ^ simple 
power-law distribution gives a reasonably good description of 
the data: 

(Rcloud) - RLuA ' With ^Sd < ^cloud < RfCd , (7) 

which implies a mass distribution: 

^ (Mcioud) - mI^J^'' , with < A^cioud < M^^^ (8) 

where Mdoud is the mass of a single cloud. This distribu- 
tion is defined between some minimum and maximum sizes. 



We generate a population of clouds as follows: in each rest- 
equivalent width bin we randomly picked the sizes of A'cioud 
clouds drawn from and calculated their ionization struc- 
ture and their H I column densities. The total H I column 
of all clouds was then summed taking into account the scat- 
ter caused by intercepting a spherical cloud at various posi- 
tions. This procedure resulted in one simulated point in Fig. 
[3] which was then repeated covering the entire rest-equivalent 
width range spanned by the data. As we are only interested 
in modeling the gas contribution originating from the halo of 
a galaxy, we do not model the high HI column density values 
with A^Hi > lO^^cm^^ and for which a disk is a likely contrib- 
utor (Zwaan et al. 2005). We note that our estimates for the 
total mass of cool gas in the halos of galaxies will change by 
a factor of order unity by including the upper end of the col- 
umn density distribution in our model. Given our density es- 
timates from photoionization calculations and the constraints 
on the mean size of clouds then a good fit for the Rao et al. 
dataset is obtained for ^Xlj = 0.3 kpc, R'^H^ = 1.5 kpc and 
j3 = —2.5 (corresponding to /3 = —0.5 for the observed dis- 
tribution). Our approach allows us to reproduce the main ob- 
servational trends summarized in Fig. |4] 

The gas temperature in our model is 10^ K and so the gas 
pressure P/kg ^ 200 cm^^ K (i.e., similar to the value used by 
e.g.. Mo & Miralda-Escude 1996 and Gnat & Sternberg 2004 
but somewhat lower than that implied by Fukugita & Pee- 
bles 2006). The implied cloud mass distribution, 0(Mcioud) °^ 
^doud' with M-Zi = 6 X 10* Mo and = 7 x 10^ Mq. 

Clouds are Jeans stable provided that the total mass (dark 
and baryonic matter) is Mdoud < 10^ Mq. The clouds are 
also relatively stable with respect to hydrodynamic instabil- 
ities (such as Klevin-Helmholtz) and evaporation (Mailer & 
Bullock 2004). 

It is interesting to note that our typical Mg II absorbing 
cloud hydrogen is mostly ionized. Most of the metals are 
in singly or doubly ionized configuration. Nevertheless, the 
large total gas columns result in higher ionization levels hav- 
ing a non-negligible column density. Specifically, the column 
density of C IV in our model can exceed lO'^ cm^^ which 
translates to an optical depth of order unity and is therefore 
consistent with a picture in which C IV and Mg II absorp- 
tion have a common origin. The uniform density, cool clouds 
considered here cannot account for the strong intervening 
O VI A 1035 systems (see Danforth & Shull 2005, Bergeron 
& Herbert-Fort 2005) and an additional component/physics is 
required to explain them (e.g.. Mo & Miralda-Escude 1996, 
McKee & Begelman 1990 and references therein). 

In our analysis we have assumed that all clouds are exposed 
to a similar flux level regardless of their location in the halo. 
Nevertheless, the effect of shielding of the ionizing flux by 
other clouds could be, in principal, important. To estimate 
the magnitude of this effect let us consider a cloud located 
50 kpc from the galaxy. Such a cloud would see an effective 
column density of a few x lO'^ cm^^ (corresponding to a few 
clouds). When shielding is included as well as the contribu- 
tion from the central galaxy, we find an overall decrease in 
the photoionization rate by a factor ^ 2 compared to the non- 
shielded case. 



2.4. Summary & limitations of the model 

As stated above, our model is aimed at providing a frame- 
work for studying the effects of a quasar radiation field on the 
absorption and emission properties of the corresponding sys- 
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Fig. 5. — The extended cool gas component considered in this work within 
the framework of the massive corona model by Fukugita & Peebles (2006; see 
their paper for the assumptions and physical motivation behind each compo- 
nent). The total mass in cool gas is comparable to that of the hot component 
within the virial radius. Given the condensed form of the extended cool com- 
ponent it occupies only a small fraction (~ 1%) of the entire volume. 




Fig. 6. — The ionization fraction of the abundant elements as a function 
of the ionization parameter. U . In this case the radiation field of the quasar 
overwhelms all other sources within ~ lOA/pc of the quasar. Gas with the 
same density and distance as a typical Mg II absorbing cloud would be highly 
ionized. 



TABLE 1 
Model summary 



Cloud parameters 


Gas mass 


6 X 10-* < MdoLid < V X IO^Mq 


Radius 


0.3 <i?doi,d < 1-5 kpc 


Density 


Pcloud = 02 atom cm^^ 


Metallicity 


0.1 Zg, 


Halo parameters (20 kpc < r < ry) 


Cool gas density profile 


p « r^^ 


Dark matter velocity dispersion 


a = 140 km/s 


Virial radius 


r„ = 282 kpc 


Virial mass 


M ~ IO'^Mq 


Cool gas Mass" 


Mcool ~ 10"Mq 


Number of clouds" 


~ 7 X 10"* 



"for a covering factor of unity and within r,. (see section 2.4) 



tem. We have calibrated it using observational constraints on 
the spatial distribution of Mg II clouds within the halo of ^ L* 
galaxies and the corresponding HI column densities, and us- 
ing photoionization calculations we have shown that the im- 
plied cloud parameters are consistent with the existence of 
Mg II gas. We summarize the main model parameters in Ta- 
bleffl 

This model leads to a picture in which L* halos are filled 
with gaseous clouds which are partially ionized and suffi- 
ciently optically thick to be detected in both low- and high- 
ionization (e.g., C IV) absorption lines. These clouds have 
sizes of order 1 kpc, masses of ^ lO^Mp^ and particle densi- 
ties ~ 10^^ cm^^, amounting to a cool gas mass of ^ IO^'M© 
within the virial radius, i.e. ^ 280 kpc. 

The results presented above assume that cool gas clouds are 
distributed throughout the entire volume of the halo up to the 
virial radius. Nevertheless, the virial radius has no clear phys- 
ical connection to the scale, rcooi, up to which cool clouds 
exist. In particular, the value of rcooi is expected to depend 
on the (unknown) origin of the cool gas. For example, if the 
gas is ejected from the galaxy rcooi is likely to be smaller than 
the case in which cool gas is accreted from the inter-galactic 



medium. Moreoever, different values may be obtained if the 
cool gas condenses out of the virialized halo gas (e.g.. Mo & 
Miralda-Escude 1996). At present, little is known with confi- 
dence about the origin of such cool gas and observations pro- 
vide only very loose constraints, hence r^ooi is ill-determined. 
That being said, we find that the mass of cool gas in galaxy 
halos depend relatively little on rcooi for the relevant scale 
range. Specifically, to reproduce the covering factor statis- 
tics on large scales, we obtain Mcool 3 (4) x lO'^M© for 
'"cool = 50 (100) kpc. Clearly, the mass of cool gas in galaxy 
halos is considerable for all plausible values of rcooi and a cov- 
ering factor of order unity. In what follows we therefore as- 
sume the gas occupies the entire virial volume (i.e., rcooi = '"vO- 

We note that, given our model assumptions and calibration 
constraints, the overall mass of cool gas, M, scales with the 
mean cloud size, (/?cioud) and is inversely proportional to the 
assumed gas density. It must be emphasized that the total 
mass of cool gas depends on the assumed covering for absorp- 
tion, i.e. a quantity still debated in the literature. In the present 
study we have used a covering factor of unity for strong Mgll 
absorption within ^ 50 kpc, as suggested by Steidel et al. 
(199 7). Ho wever, other studies claim lower values (~ 0.2 
from iBech told & Ellingson ( 1992) and ^ 0.5 from Tripp & 
Bowen 2005). As the total mass of cool gas scales linearly 
with this quantity, we might be overestimating the mass if the 
actual covering factor is not unity. 

According to our model, the mass of cool gas in galaxy ha- 
los scales linearly with halo size. Specifically, the mass within 
30kpc is Mcool < IO'^Mq and, by extrapolation to smaller 
scales, is comparable to that estimated by Martin (2006) for 
ultra-luminous infrared galaxy outflows. In Fig. I2.4l we show 
the different components of the average mass profile of late- 
type galaxies proposed by Fukugita & Peebles (2006), and 
show that including the distribution of cool gas proposed in 
the present paper is only a modest departure from it. 

So far we have concentrated on modeling the gas within the 
halo of an L* galaxy with My,> = 2.5 x 10 M©. It is interest- 
ing to note that our model may be scaled to different galaxy 
luminosities and masses. Indeed, analyses by Steidel et al. 
(1997) and Guillemin & Bergeron (1997) have indicated that 
the normalization of the Mg II rest equivalent width follows a 
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Fig. 7. — The thermal and ionization properties of halo gas {/m = 5) illuminated by a quasar. Left: The temperature profile for a constant density, optically 
thin medium with density of 0.02 cm"' (see text) for several values of quasar luminosity. Clearly, the temperature drops with distance as the ionization parameter 
decreases {<x r^^). The radiative effect of Seyfert galaxies on their environment is more localized than that of bright quasai's. Right: The ionization structure as 
a function of distance from a L^f, = 1 quasar As expected, the ionization level decreases with distance due to the dilution of the radiation field. The fraction of 
Mg II recovers only on Mpc scales while closer to the quasar (on ~ lOOkpc scales) C IV is expected to show absorption signatures while on smaller scales still 
20kpc) O VI would show up in the spectrum. The gas is very highly ionized within ~ 30kpc and could be detected in the X-rays. For quasars with different 

luminosities the abscissa should be rescaled by L45 . Thus, brighter quasars would ionize their environment to lai'ger distances. The dashed line con'esponds to 
the ionization fraction of H I for the case of L46 = 0.1. The hump beyond ~ 50kpc is the result of self-shielding (see text). 



relation 



Wo(r,L)=Wo(r)x 



(9) 



with Y — 0.1 — 0.4. Such a scaling has been reported in the 
luminosity range 0.3 — 5L* and at redshifts 0.5 < z < 1. Ex- 
trapolating it further is therefore uncertain. The above relation 
may be understood within the framework of our model as a 
dependence of the number of clouds, hence mass of cool gas, 
on luminosity. Denoting the mass of cool gas in L* galaxies 
by M* then the mass of cool gas in halos of galaxies with lu- 
minosity L may be written as M{L) = fM{L)M* where fm is 
a dimensionless scaling factor. 

3. QUASAR HALOES 

We now investigate the effects of a quasar on the gas distri- 
bution within a halo and discuss the observational signature in 
absorption and emission. Our main assumption is that quasar 
halos are analogous to those of non-active galaxies. As indi- 
cated by recent studies (e.g., Jahnke et al. 2004), the hosts of 
bright quasars seem to correspond to a fewxL* galaxies. Fur- 
thermore, a recent study by Serber et al. (2006) seems to indi- 
cate that bright quasars reside in over-dense regions implying 
that their dark matter halos are several times those typical of 
L* galaxies. Motivated by these studies we shall allow /m to 
deviate from unity. We show below that many observational 
constraints may be accommodated by using /m ^ 5 — 10 for 
the mass of cool gas in quasar halos. 

3.1. Thermal and Ionization Structure 

In what follows we compute the time-independent ion- 
ization and thermal structure of a cloud exposed to a con- 
stant quasar ionizing flux. In the Appendix we present time- 
liepeniienf calculations and show that the steady-state assump- 
tion suffices for the luminosity range and the relevant spectral 
features considered here. 

In modeling the spectral energy distribution (SED) of the 
quasar we have used the continuum defined by Sazonov, Os- 
triker, & Sunyaev (2004) which is shown in Fig. |2] Ionizing 
sources other than the quasar can be safely neglected. We start 



by introducing the ionization parameter, U, which is the ratio 
of photon density to particle number density and is a measure 
of the ionization level for photoionized gas. For the chosen 
SED we obtain that 



f/~0 



6L46 



-1 



50kpc 



(10) 



where the quasar luminosity is L = 10'*^L46 erg s"'. The 
ionization structure of several commonly detected ions as a 
function of U is shown in Fig. |6] By combining these re- 
sults with the density estimates from our model (§2), we can 
compute the ionization level and temperature as a function 
of the radius from the center of the galaxy. The results for 
the temperature profile are shown in Fig. [7] As expected, 
the quasar has a substantial effect on the thermal state of the 
inner parts of the halo and can heat the (initially cool) gas 
to high temperatures by means of photo-absorption. The ex- 
tent to which the Mg II halo is ionized depends mainly on the 
quasar luminosity. For example, a quasar with L46 = 1 heats 
gas within the inner 20kpc to some 10^ K with gas at 50kpc 
being ionized to a few xlO^K. Objects with Seyfert-like lu- 
minosities (L4(5 < 0.1) have a smaller effect on their immedi- 
ate environments and heat the gas within their inner 20 kpc to 
a few X 10^ K with clouds at larger radii remaining relatively 
cool. In contrast, very luminous quasars can heat their 100 kpc 
environment to temperatures of order 10^ K with gas within 
the inner 10 kpc being able to reach the Compton tempera- 
ture. The effect of a luminous quasar is not necessarily lim- 
ited to halo gas and low density gas may be affected on much 
larger (Mpc) scales. Our calculations indicate (Fig. ^ that 
for clouds to remain relatively cool within the central 30 kpc 
of a L46 = 1 source, their density should be ^ 10^ cm"^, i.e. 
considerably higher than the values considered here and those 
thought to be representative of Mg II clouds (c.f., Stockton et 
al. 2002). 

The ionization structure closely follows the temperature 
profile. Fig. [T] shows the ionization fraction of several ions 
as a function of radius in the halo for the case of L46 = 1. 
Clearly, gas within the inner 100 kpc is devoid of Mg II, with 
C IV O VI, and O VIII being most abundant; the ionization 
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Fig. 8. — The geometry considered here whereby a quasar emits ion- 
izing radiation into a cone with opening angle 9. Cool gas clouds which 
are exposed to the ionizing radiation (circles) are heated while those that 
are shielded from the ionizing radiation remain cool (filled circles) and their 
properties resemble those of halos associated with non-active galaxies. Off- 
center lines of sight will intercept cool gas in the transverse direction, while 
those that point radially toward the quasar will mostly see ionized gas (see 
text). 





0, 
I 

go. 

03 

o " 






f„=1 ; 8=45" 

f =1;e=70° 
M 

f. =5; 9=45° 
f,,=10;e=70° 

M 




50 100 150 200 250 

kpc 

Fig. 9. — The covering factor for strong Mg II absorption as a function of 
the impact parameter for lines-of-sight passing through the halo of a quasar 
(see figure[8). Thick (thin) lines coirespond to the covering factor for systems 
with Wo > lA (Wo > 0.6A). Four halo models are considered that differ by 
the mass of cool gas within the virial radius (denoted by fu; see text) and 
the opening angle, 9, of the quasar radiation field. As expected, the covering 
factor declines with increasing impact parameter due to the decreasing cloud 
density with distance. The covering factor for strong absorption is sensitive 
to both 9 and . As 9 increases, more gas is ionized and the halo becomes 
transparent to Mg II absorption. As /m rises, more clouds are intercepted 
along lines-of-sight through the halo. As shown, different combinations of 
fm and 9 can result in a similar covering factor for strong Mg II absorption 
(see text). 



fraction of the latter peaking on the smallest scales. For a 
/m — 5 halo with a L46 — 1 quasar (e.g., Jahnke et al. 2004 
and equation |9]l, the gas remains optically thin above the 
Lyman edge. For different source luminosities the abscissa 
should be rescaled by L^^. For low quasar luminosities 
(L46 < 0.1) even the inner regions of the halo remain at low 
temperatures and the gas is less affected by the quasar on 
~ 100 kpc scales (see Fig.|7ll. 

The active galactic nuclei (AGN) unification scheme (An- 
tonucci 1993) suggests that quasars emit light in a double 
cone geometry with an opening angle (see Fig. [8]l. Such 



an anisotropy in the radiation field will be imprinted in the 
properties of the gas distribution around the quasar. In the 
next sections we discuss the observational implications, both 
in absorption and emission. 

3.2. Quasar-halo absorption properties 

Clouds of gas that lie within the ionization cone of the 
quasar would be heated and ionized, and would follow the 
ionization structure presented in Fig. [T] Therefore, the line- 
of-sight to a quasar would be largely devoid of associated 
Mg II absorption due to the halo clouds modelled in this pa- 
per Nevertheless, such highly ionized clouds are expected to 
give rise to associated absorption features from species like 
H I, C IV, and O VI - depending on their ionization structure 
which is determined by the quasar luminosity and their dis- 
tance from the source. In contrast, lines-of-sight probing the 
quasar halo at finite impact parameters from the quasar would 
intercept cool gas which is not affected by the quasar radiation 
field and is therefore likely to show absorption features due to 
low ionization species as those seen around non-active galax- 
ies. The number of Mg II clouds intercepted by a given sight- 
line depends not only on the total number of clouds within the 
halo (which, for non-active galaxies, depends on the galaxy 
luminosity; see equation 9) but also on the opening angle of 
the quasar. 

Given the above model, we now discuss the observable 
properties of quasar haloes in terms of their expected ab- 
sorption signatures. In a recent work we have used pairs of 
quasars at different redshifts and investigated the presence of 
Mg II absorption around the foreground objects (Bowen et al. 
2006). In 4/4 cases we have detected Mg II absorption with 
Wo > O.5A in the spectrum of the background quasars with a 
redshift matching that of the foreground quasar This suggests 
that strong Mg II absorbers are common in the transverse di- 
rection of quasars. We have extended our analysis to a larger 
sample of systems and our new results (Bowen et al. 2007) 
indicate that the presence of strong Mg II absorption around 
quasars with a covering factor close to unity within 100 kpc, 
with a rapid decline on larger scales. Such behavior is sim- 
ilar to that of Mg II gas which is detected around galaxies. 
Interestingly, the high covering factor for Mg II gas in the 
transverse direction seem to be in contrast with smaller val- 
ues (-- 10%; Aldcroft et al. 1994) derived for Mg II absorp- 
tion along the line-of-sight to these objects (see Fig. [8]). As 
we show below, this surprising finding is naturally explained 
by our model. 

Quasar radiation is thought to be emitted into cones of 
opening angle 9 (Antonucci 1993). As such, is a measure 
of the mass of cool gas within the halo which is exposed to the 
ionizing field and determines the mass of the remaining cool 
gas. The connection between 9 and Wo is therefore apparent 
within the framework of the model (see figure[8]l. Model pre- 
dictions for the absorption covering factor as a function of the 
rest equivalent width, the impact parameter, and the mass of 
cool gas, /m, are given in Fig. |9l These results are analogous 
to those of galaxies but for the additional dependence on 9. 
As shown, quasars that emit more isotropically have less cool 
gas at any impact parameter due to ionization effects. More 
massive envelopes of cool gas have a larger number of clouds 
and so result in broader Mg II troughs. Clearly, different com- 
binations of /m and 9 can yield similar covering fractions. 
At present, neither the mass of cool gas in quasar halos nor 
the opening angle of quasar ionization cones are well deter- 
mined with recent surveys indicating that of 9 may depend on 
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Fig. 10. — The predicted covering factor for associated absorption where 
lines-of-sight probe the ionized part of the halo. Covering factors for several 
high ionization Unes are considered correspnding to: Hi A1216 (dash-dotted 
line), IVIg II A2800 (dotted line), C IV A 1548 (solid line), and O VI A 1035 
(dashed line). The results are shown as a function of the quasar luminosity 
for L46 = 0.1 — 10 (columns) and mass in cold gas with = 1 — 10. As 
shown, the covering factor is more sensitive to the luminosity rather than the 
mass in cold gas. The covering factor follows the exponential cutoff of the 
Poisson distribution for the respective lines. Low luminosity quasars have a 
higher covering factor for low ionization lines since the halo is less ionized. 
For bright objects the halo is more ionized and low ionization hnes become 
weaker. 

the quasar luminosity and is in the range 45° < < 70° for 
bright quasars (Willott et al. 2000, Treister & Urry 2005). Fig. 
|9] shows that /m — 5 — 10 is required to explain the high oc- 
currence of transverse strong Mg II absorption around quasars 
within 100 kpc in the Bowen et al. (2006, 2007) sample. 

We can further compare the predictions of our model to ob- 
servations by considering the fraction of quasars showing as- 
sociated absorption, i.e. absorption arising along the quasar 
line-of-sight close to the quasar redshift (Fig. O. As men- 
tioned above, in such a case, low-ionization species are not 
expected to be abundant due to rapid ionizations by the quasar 
radiation field. However, the same clouds are expected to 
give rise to higher ionization absorption lines. In Fig. [TOlwe 
present predictions for the covering factor of H I, Mg II, C IV, 
and O VI lines as a function of halo mass and quasar quasar lu- 
minosity. The covering factor for all lines is seen to decrease 
exponentially with increasing Wq due to Poisson statistics. As 
shown, the covering factor for absorption is more sensitive 
to the quasar luminosity than to the mass of cool gas in the 
halo. In particular, the halos of more luminous quasars are 
more ionized resulting in a lower covering factor for strong 
low ionization systems. As expected, more massive halos pro- 
duce stronger absorption due to the larger number of clouds 
intercepted along the line-of-sight. Predictions appear to be in 
qualitative agreement with recent surveys indicating a cover- 
ing factor of 30% for strong C IV 1548, 1550 absorption 
(Vestergaard et al. 2003). Our model predicts little associated 
Mg II absorption due to halo gas. We note however, that our 
model does not include the, possibly important, contribution 
from absorption by gas in the immediate vicinity of the black 
hole (sometimes termed intrinsic absorption). That said, our 
model is consistent with the notion that most low ionization 
associated systems are intrinsic to the quasar. Additional data 
concerning the covering factor of various ionization stages as 
a function of quasar luminosity (and perhaps host mass) are 
required to test our model. 

We note that a complementary configuration to the one used 
in the Bowen et al. (2006, 2007) survey is provided by radio 



galaxies. For these objects our line-of-sight passes outside 
the quasar radiation cone. Measurements of large-scale H I 
absorption by Van Ojik et al. (1997) have shown that typical 
column densities through the neutral part of the halo are of 
order < 10^^ cm^^. Such values are in qualitative agreement 
with our model and with the Rao et al. (2006) dataset (see 
§2). A more detailed comparison is unwarranted at this stage. 

3.3. Quasar-halo emission properties 

As the gas is photo-excited by the radiation field of the 
quasar, clouds would not only absorb but also emit and scatter 
continuum and line photons. In this section we calculate the 
expected emission line flux from the halo. 

Line emissivity depends on the ionization state of the gas as 
well as on the optical depth through which line photons need 
to propagate, escape, and reach the observer In calculating 
the emissivity of lines we have assumed that once line pho- 
tons escape individual clouds, they reach the observer; i.e., 
we assume no photon scattering or absorption among differ- 
ent clouds. This approximation is adequate if clouds along 
our line-of-sight have little overlap in velocity space, which is 
the case considered in this paper The fraction of line photons 
escaping from a single cloud was calculated using CLOUDY 
using the escape probability method for a static cloud. The 
resulting emissivity per unit mass as a function of distance 
from the ionizing source is shown in Fig. [TT|for Lya emis- 
sion line. As shown, averaged over volume, Lya emissivity 
is relatively insensitive to the luminosity of the central object. 
This results from the gas temperature having only a weak de- 
pendence on luminosity for the relevant parameter range and 
from hydrogen being fully ionized. Emissivity is suppressed 
on small scales near very luminous sources due to suppression 
of the recombination rate at high (^ 10^ K) temperatures. 

We have also studied the emissivity of other lines and, in 
pai-ticular, that of [O III] A5007 forbidden line (Fig. [TTT ). 
The line emissivity in this case is very different from that 
of Lya and other hydrogen-like lines since the recombina- 
tion rate is much lower (O III and O IV abundances in the gas 
are negligible). Relatively neutral gas, which is required for 
efficient [O III] A5007 emission, becomes abundant only on 
large scales and even then only for relatively low luminosity 
quasars. As quasars get brighter the peak O III emissivity is 
shifted to larger scales where cooler gas is found. We em- 
phasize that, unlike Lya emission, [O III] A5007 emission is 
sensitive to gas composition and is roughly linear with it. 

Our calculations allow us to predict how quasar halos would 
appear for different emission lines. We show an example for a 
Lya image in Fig. [TT]for the case of = 70°, /m — 5, L46 — 1 
and assume that our line-of-sight is along the symmetry axis 
of the quasar. The image is symmetric with the surface bright- 
ness profile being roughly ' from the center (neglecting 
the effects of scattering; see below). We note that asymmetric 
images are obtained if the axis of symmetry is at an angle to 
our line-of-sight, as likely to be the case in general. 

Different lines have different emissivity profiles and hence 
different luminosities and inferred nebular sizes. Table 2 lists 
a few of the more important lines which could, in principal, 
be observed in quasar haloes. Lya, C IV A 1548, and O VI 
A 1035 are the strongest transitions with the latter two lines 
depending roughly linearly on the Metallicity. The typical 
nebular extent (from which most of the flux is emitted) of low 
ionization lines is > 100 kpc. This is in contrast to the appar- 
ent size of high ionization lines which come from compact, 
highly ionized regions of the halo (e.g., the nebular extent of 
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TABLE 2 

Emission Line Properties of Quasar Nebulae 



Line ID 


L 


Size 


^scal.miix 




[erg s"'] 


[kpc] 


[] 


HiLya 


3 X 10"^' 


150 


0.25 


VI A 1035 


4 X lO''^ 


80 


0.25 


NV A1238 


2 X 10« 


150 


0.25 


CIV A 1548 


2 X 10^2 


200 


0.5 


HI Ha 


1 X 10« 


100 




HlH/5 


8 X lO''^ 


150 




C III Mil 


7x 10"' 


200 




O VIII A 18.97 


4x 10"' 


20 


7 


[OIII] A 5007 


2 X 10"' 


200 





Luminosities for doublets are the sum for individual transitions. Size is defined 
as the extent of the emitting nebula which contains 90% of the flux for a given 
transition. The last column is an order of magnitude estimate for the maximum 
luminosity in scattered fight (see text). Transitions that do not scatter are marked 
with "-" while those for which the broad emission line flux is very uncertain 
are marked with "?" (see text). L46 = 1,6 = 70° and /m = 5 are assumed. 



O VIII X-ray lines is > lOkpc for L46 = 1). 

3.3.1. Scattering 

In addition to line emission there is scattering (i.e. ab- 
sorption followed by emission) of the photons emitted by the 
quasar. Here broad emission line photons can be scattered off 
gas clouds in the halo provided the neutral hydrogen fraction 
is large enough. Owing to the low gas density, the population 
of excited levels is negligible and so only resonance lines are 
able to scatter efficiently. To correctly calculate the surface 
brightness profile due to photon scattering one requires de- 
tailed knowledge of the 3D distribution of clouds in real and 
velocity space. In addition, the fraction of photons reaching 
the halo depends on the amount of intrinsic absorption and the 
relative velocities of the emission line and halo gas. Detailed 
calculations should include the effect of multiple scatterings 
even in this relatively optically thin medium. 

Despite the overall complexity of scattering problems, it 
is relatively straightforward to estimate the maximum scat- 
tered flux expected for different lines. This is simply the 
flux in the broad emission line with velocity dispersion, C7/ 
which is intercepted by halo clouds with velocity dispersion, 
(7c given by roughly (Jc/<7i — 10% (here we used (7/ ~ 2000; 
e.g., Vanden Berk et al. 2001). Given our quasar SED and 
taking the mean rest equivalent width for Lya to be ^ 90A 
(e.g., Zheng et al. 1997; Vanden Berk et al. 2001) we find 
that the the maximum scattered Lya luminosity, 1!^^]^™ ~ 

2 X 10"*^ [1 - sin(0)]L46 erg s"^ Table 2 lists the results for 
other lines based on the same line-of-reasoning and shows that 
the contribution of scattered broad line region photons to the 
large scale nebular emission is probably less than 50% for all 
prominent lines (and forL46 < 10^). 

To better evaluate the flux due to broad emission line pho- 
tons scattering off halo gas, we define an effective optical 
depth over the (dynamically broadened) absorption line pro- 
file of halo gas, Tes{r) — C(r)T(r) where C(r) is the geometric 
covering factor of clouds at a given distance (i.e., the number 
of clouds intercepted at some distance interval) and T(r) is 
the line optical depth of a single cloud at position r. In our 
model, C(r) is large on small scales due to the clouds den- 
sity peaking toward the center yet, T(r), is low since the gas 
is more ionized. For Lya, Tgff is decreasing with decreasing 
r. On scales larger than ^ 30kpc, clouds become opaque to 
Lya and Teff °^ C °^ r^^. We have approximated the problem 



by running monte-carlo simulations (e.g., Zheng & Miralda- 
Escude 2002; Z. Zheng, private communication) of photons 
emitted by a point source and being scattered off halo gas 
with Tgff being of order a few and having the radial depen- 
dence described above. We find that the surface brightness 
of the scattered photons is 0= b^^-^ and therefore that most 
of the scattering occurs in the inner nebular regions (within 
~ 30kpc). The results are shown in the inset of Fig. [TTIfor 
the case of L46 — 1 . We note, however, that these provide only 
order-of-magnitude estimates since the gas kinematics in the 
central regions of the halo is rather uncertain. 

3.3.2. Detectability 

It is interesting to examine the detectability of line emitting 
nebulae around quasars within the framework of our model. 
Here we focus on Lya nebulosity and study its observed 
luminosity and size as a function of the limiting flux of the 
observation (see Fig. [T2]i. As before, fu = 10, = 70° and 
L46 = 10 are assumed. We also examine these quantities 
for several redshifts assuming the observed nebulae at 
those redshifts are similar so that the surface brightness is 
°^ (1 Our calculations indicate that the limiting flux 

required to observe extended Lya emission around z = 1 
quasars is better than ^ 10^'^ erg s^' cm^^ asec^^. In 
particular, our model suggests that the observed nebular 
extent is sensitive to the limiting flux with a factor two 
increase in sensitivity allowing to detect emission from an 
order of magnitude larger volume. At low-z, a limiting flux 
of a fewxlO^'** ergs^' cm^^ asec^^ is required to probe 
the full extent of the nebulosity and correctly estimate its 
luminosity. Once the appropriate sensitivity is reached, our 
model suggests that large scale emission line nebulae in 
general, and Lya nebulae in particular, should be a rather 
common phenomenon around quasars. More massive halos 
or quasars having a wider ionization cones are better emitters 
and their halos can be observed to higher-z. Our calculations 
suggest that Lya scattering can have a substantial contribu- 
tion to the observed luminosity (compare models B,C and D 
in figure [T2b . As expected, halos with lower mass of cool 
gas {Jm = 1) are poor emitters (compare models A and B in 
figure [12]). Similar considerations also apply for metal lines 
whose detection can be used to put more stringent constraints 
on the ionization level of the emitting gas as well as estimate 
its mean metallicity. For example, the N V A 1238 emission 
line could have a flux comparable to that of Lya for >solar 
metalHcity gas (see table 2). 

3.3.3. Comparison with observations 

Extended narrow-line emission around quasars and radio 
galaxies has been observed for more than two decades and 
the accumulation of data has shown that giant Lya nebulae 
can be found around QSOs, up to scales of ^ 100 kpc (e.g., 
Wampler et al. 1975, Stockton 1976, Bremer et al. 1992, 
Boroson, Oke, & Green 1982, Bergeron et al. 1983, Boro- 
son, Persson, & Oke 1985, Stockton & MacKenty 1987, Reu- 
land et al. 2003, Labiano et al. 2005). Early studies focused 
mostly on radio loud quasars but more recent ones showed 
that radio-quiet quasars also show extended Lya emission on 
large scales. It is interesting to note that Hu et al. (1991) 
observed a sample of 7 radio-quiet quasars down to a limit- 
ing flux of 2 X 10^'^ergs^'cm^^ and found no extended Lya 
emission. However, subsequent and more sensitive obser- 
vations showed that quasar nebulae are commonly observed 
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Fig. 11. — The Lya (solid) and [O III] A5()07 (dashed) emissivities for various color coded luminosities. While the Lya emissivities of the nebulae at large 
radii are similar, the emissivity decreases with increasing quasar luminosity on smaller, < 30 kpc scales due to ionization effects. The [O III] A5007 emissivities 
are lower than those of La and have their maxima at larger radii where cool gas is more abundant. Right: The Lya halo surrounding an L^f, = 1 (fu = 5,9 = 70°) 
quasar at z = 1 with colors reflecting the flux level. Also shown in an inset is the surface brightness profile including (solid line) and excluding (dashed line) 
scattering of Lya photons from the broad emission line region of the quasar (see text). 



around radio-quiet quasars (Weidinger et al. 2005, Bremer 
et al. 1992., Steidel et al. 1991, Hu et al. 1996, Petitjean 
et al. 1996). These authors report Lya luminosities in the 
range 7 — 58 x 10'*^ergs^^, values which are largely consis- 
tent with our predictions as shown in Fig. 12. Christensen 
et al. (2006) report Lya luminosities from z ^ 3 nebulae 
which are of order ^ 10"*^ erg s^^ and extend up to a me- 
dian value of 15 kpc (with some extending up to 60 kpc) for 
a flux limit of ^ 10^'^ erg cm^^ s^' asec^^. This is consis- 
tent with quasar halos having a considerable amount of gas 
in cool form (model A in figure [12] corresponding to /m = 1 
with no line scattering is excluded by the data; c.f. Binette et 
al. 2006). Christensen et al. place limits on C IV line emis- 
sion down to fluxes 10 times lower than Lya. In our model 
this is roughly the predicted flux ratio of the lines for 0. 1 solar 
metallicity. More detailed modeling is, however, beyond the 
scope of this paper 

The statistical properties of quasar halos and their under- 
lying physics are poorly known and few models have been 
devised to account for them (e.g., Heckman et al. 1991). 
Common explanations for the existence of such line emission 
nebulae include gas accretion onto the central object (Haiman 
& Rees 2001) and matter ejection by the quasar and its host 
galaxy (e.g., Heckman et al. 1991). Moreover, the interpre- 
tation of the data is controversial: while several teams have 
reported high pressures of order P ~ 10^^'cm^^ K for radio 
loud quasars (RLQ), a study by Fu & Stockton (2006) points 
to lower pressures, of order lO^cm^^ K. These estimates 
come from compact emission regions of [Oil] and [Olll] 
lines and from absorption by excited states (e.g., Hamann et 
al. 2001). In our model for the halo of an L* galaxy the 
pressure is even lower, P ~a fewxlO^ cm^^. It is possible 
that RLQ are different from radio-quiet quasars (RQQ) given 
their overaU higher Lya luminosities (e.g., Christensen et al. 
2006) and their apparently large X-ray and [OIII] emission on 
small scales (e.g., Crawford & Fabian 1989). Our model pre- 
dicts relatively faint and diffuse [O III] emission from the halo 
which seems to be at odds with recent detections of [O II] and 
[O III] emitting filaments around RLQs. Deep optical obser- 
vations of low-z RQQs may help to determine whether there 
is a relatively dense and filamentary structure in quasar halos 
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observed nebular luminosity [erg s ] 



Fig. 12. — The observed Lya nebula size (containing 90% of the to- 
tal obsrved flux) vs. the observed luminosity for a halo illuminated by a 
L46 = 10 quasar with an opening angle 9 = 70°. The following models are 
considered: (A) fm = ^ with no line scattering. (B) fu = 10 with no line 
scattering. (C) /m = 10 with scattered Lya luminosity, Lf^^ = O.SL^y*'"'"'. 

(D) fM = 10 and Ll^q ~ ^Lya'""- assume no evolution in the nebular 
properties with redshift. Placing the same nebula at different redshifts (de- 
noted by different symbols; see legend) decreases its observed luminosity 
and size for a survey with a given limiting flux (denoted by colors in units 
of erg s^'cm^^asec^^). Clearly, unless the nebula is close and/or the sen- 
sitivity is high, only the brighter inner part of the nebula is observed which 
results in an under-estimation of its luminosity and size (see text). In particu- 
lar, sensitivities of order 10^'^ erg s^'cm^-asec^^ are required to probe the 
full extent of quasar nebulae at j < 2 (blue square and circle). 

in addition to the more diffuse component considered here. 

4. DISCUSSION 
4.1. Galaxies 

Our model implies considerable mass in cool gas out to 
the virial radius of L* galaxy halos at z ^ 1. As explained 
above, this assumes a unit covering factor for strong absorp- 
tion within 50 kpc of the galaxy. The lower covering fac- 
tors suggested by Bechtold & Ellingson (1992) and Tripp & 
Bowen (2005) would actually lead to an agreement between 
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our model predictions and the results of Mailer & Bullock 
(2004) who estimated the mass of cool gas in the Galaxy halo 
to be 2 X lO'^ Mq and argued for consistency with observa- 
tions of high-velocity clouds. At present, radio surveys do 
not possess the required sensitivity to probe low column den- 
sity, ~ lOOkpc diffuse gas in 21 cm emission, as predicted by 
our model. We caution, however, that it is not clear whether 
our model predictions for z ~ 1 objects are directly applica- 
ble to present day galaxies. In particular, Nestor et al. (2005) 
find evidence for evolution of cool gas properties with cosmic 
time whereby the occurrance of strong systems appreciably 
declines from z = 1 to z = 0.5. Understanding the distribution 
of cool gas around present day galaxies will require extensive 
UV spectroscopic surveys. 

The notion of massive galaxy halos has been promoted by 
Fukugita & Peebles (2006) to provide a solution to the miss- 
ing baryon problem. In their model, the gas mass up to the 
virial radius is a fewxlO'^ M©. Their total gas mass es- 
timates provide a natural explanation to the missing baryon 
problem and are consistent with a large range of observational 
constraints which are independent of those considered here. 
Their model did not consider the possible existence of cool 
gas beyond the extent of the disk. Here we have shown that 
converting a fraction of the hot plasma into cool clumps al- 
lows us to reproduce a number of observational constraints 
from absorption line studies, while not significantly changing 
the total amount of gas within the virial radius (see Fig. 5). 

More work is needed to establish how much mass there is in 
the halos of present day galaxies. On the observational side, 
better estimates of the covering fraction of cool gas around 
galaxies are required, possibly as a function of galaxy type. 
On the theoretical side, numerical simulations of galaxy for- 
mation need to include realistic cooling rates and reach high 
resolutions to reliably trace cool and condensed forms of gas. 
Understanding the properties of the dilute and eUusive, yet 
likely very massive, cool gas component in galaxy halos is es- 
sential for understanding galaxy formation and evolution from 
the early Universe, to present times and beyond. 

4.2. Quasars 

Quasars are known to reside in galaxies yet their large scale 
environments are poorly understood. The recent discovery of 
Mg II absorbing gas around those objects (Bowen et al. 2006, 
2007) suggests an interesting hnk between galaxy and quasar 
environments. This impUes that the presence of a quasar can 
be used to shed light on the distribution of matter in the halos 
of galaxies. 

We have shown that it is possible to construct a model 
which unifies several distinct astrophysical phenomena per- 
taining to intervening absorption, associated absorption, and 
quasar emission nebulae. In our model, these phenomena are 
different manifestations of the same physical entity, that be- 
ing of cool gas condensations distributed within the halo of 
galaxies. If our model is correct then the properties of cool 
gas in the vicinity of galaxies may be revealed by studying 
the influence of an active nucleus on its thermal and ioniza- 
tion structure. In particular, the study of associated absorption 
may help to constrain the density and distance of halo gas 
from the quasar. Cool gas around galaxies is notoriously dif- 
ficult to detect in emission. Nevertheless, the presence of an 
ionizing source enhances emission and renders the detection 
of circum-galactic material more accessible. Furthermore, by 
studying the nebulosities around quasars, one can obtain a 2D- 
spatialx ID- velocity picture of the gas. Mean quantities per- 



taining to the mass of the gas and its metallicity can be more 
easily deduced for individual objects. 

In addition to shedding fight on gas distribution in the ha- 
los of galaxies, it is possible to deduce important conclu- 
sions about quasar physics: the recent results by Bowen et 
al. (2006, 2007) indicate a covering factor of order unity for 
cool gas in the transverse direction around quasars in con- 
trast to the much lower covering factor for associated Mg II 
absorption. Our model implies that this is merely a conse- 
quence of the quasar uiufication scheme: gas in the radial di- 
rection to the quasar is heated to high temperatures while that 
which is at right angles to the quasar ionization cone remains 
unaffected. This is not a trivial result since, although the uni- 
fication scheme is well established for active galactic nuclei, 
it is less secure for bright quasars at high-z. Further tests of 
this conjecture may be carried out by comparing model pre- 
dictions to the occurrence of associated absorbers for a wide 
range of ionization levels as well as by studying emission line 
nebulae in detail. 

Our model can reproduce the emission and absorption prop- 
erties of quasars by requiring quasar halos to be more ex- 
tended and therefore more massive than those of L* galax- 
ies (c.f. Serber et al. 2006). Such a scaling combined with 
the radiation properties of quasars can simultaneously explain 
the anisotropic absorption properties and the Lya nebulae ob- 
served on scales reaching ~100 kpc. Finally, if the correlation 
found between the size of cool gaseous halos and galaxy lu- 
minosity (Steidel et al. 1997, Guillemin & Bergeron 1997) 
can be extended to quasars then this impUes that quasar hosts 
are a few xL*, which is in agreement with results from Jahnke 
et al. (2004). 

Better understanding of quasar emission line nebulae may 
allow us to deduce the long-term light-curve behavior of 
quasars. This will require better understanding of the mat- 
ter distribution in quasar halos which can be achieved by deep 
imaging and spatially resolved spectroscopy of quasar envi- 
ronments. The detection of orphan nebulae (where the active 
nucleus is obscure) may provide a new means for identify- 
ing type-II quasars. Combining large samples of emitting and 
absorbing gaseous nebulae may allow us to statistically con- 
strain the mass of cool gas in quasar halos as well as the open- 
ing angle of quasars. This has implications for re-ionization 
and background radiation determination. A survey of emis- 
sion line nebulae around a large number of quasars may also 
reveal environmental differences among quasar types (e.g., 
the RQQ-RLQ dichotomy). 

The detection of metal emission lines, whose presence (de- 
pending on the gas composition) is predicted by our model, 
may allow to estimate the metallicity of gas in the halo and 
test models for metal enrichment of galaxies and the inter- 
galactic/intra-cluster medium. 

5. SUMMARY 

We have presented a phenomenological model for the distri- 
bution of cool gas around L* galaxies, calibrated with a wide 
range of observational constraints from absorption fine stud- 
ies (rest equivalent width distribution of Mg II absorbers, HI 
column density, etc.). We argue that the halos of L* galax- 
ies are filled with gaseous clouds with sizes of order 1 kpc, 
masses of ~ lO^M© and particle densities ~ lO^'^cm"^. The 
total amount of cool gas within the virial radius of L* galax- 
ies depends on the covering factor for strong Mg II absorption 
(which is a matter of debate) and is likely to be in the range 
10^*^ < Mcooi < 10^^ M©, within the virial radius. 
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By assuming self-similarity, we show that the obtained so- 
lution for the distribution of cool gas around galaxies, if ap- 
propriately scaled, can reproduce the properties of cool gas 
seen quasars. Our model simultaneously provides an expla- 
nation for the La nebulosities observed observed on ~ 100 
kpc scales around quasars as well as the main properties of 
the cool gas seen in absorption around quasars. 

Comparison of model predictions with future surveys will 
shed light on the missing baryon problem and will deepen our 
understanding with respect to galaxy formation and quasar ac- 
tivation. 
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APPENDIX 

Here we describe the time-dependent ionization and thermal structure of an initially cool cloud which is exposed to a quasar's 
ionizing radiation. As the quasar phenomenon is short lived compared to the Hubble time (the quasar lifetime, Tg ^ 10^ years; 
e.g., Worseck & Wisotzki 2006), it is instructive to first consider a few relevant timescales of the problem; 

• The dynamical timescale of a halo of size ^haio at a virial temperature, Ty (corresponding to some sound speed, Vj) is 

^halo ,„9 '■ f V^''^ 

"''^"^^"^^TookF^llo^j 

• The sound crossing timescale for a cloud of size 2/?(.ioud at temperature T is 

27?cloud ,„8^cloud f T \ 

—^'^^[wk) 

• The photoionization timescale for ion X is 

where ax is the cross-section, E the energy, Le the flux per unit energy, r the distance from the quasar. We assumed a mean 
ionizing photon energy of £ = 10 eV and an ionization cross-section of 10^'^ cm^ (roughly that of Mg II). 

• The recombination timescale for ion X is 

t£,^ ^—-^5xlO^( :r] ' yr (4) 

nax{T) V0.02cm-3y ^ ^ ^ 

where n is the number density of atoms and ax the recombination coefficient (taking only the radiative term for Mg II). 

The photoionization and recombination timescales are the shortest in the problem. The cloud and halo dynamical timescales 
(approximated by the sound-crossing time) are probably longer than the quasar lifetime, t^; hence, while the ionization structure 
of the clouds can vary rapidly (see below), the structure of the halo and embedded clouds remains approximately constant during 
the quasar lifetime. We note, however, that the sound-crossing timescale for the smallest clouds in our model could be as short 
as ~ 10^ years if their temperature rises to high values once the quasar activates. We choose to ignore this complication here 
since (a) the quasar lifetime is poorly known and (b) the evaporation timescale is a few times the sound-crossing timescale (e.g., 
Bertoldi 1989). Also, due to the spectrum of the cloud's size-distribution, the evaporation of the smallest clouds will make little 
difference to our final conclusions. 

In addition to the above analytic approximations, we have numerically computed several time-dependent photoionization and 
thermal models. This is necessary to quantify just how important the effects of the varying quasar flux are compared to time 
independent calculations and to see which ions are affected most. The algorithm used for these calculations is fully described in 
Gnat & Sternberg (2006). 

At the beginning of the calculation we assume that the cloud is exposed to the meta-galactic field and is in ionization and 
thermal equilibrium. We then turn on a quasar and assume a step function description for its light curve, i.e., the quasar if off for 
f < and on at f > 0. The quasar is assumed to maintain a constant flux level to f ^ oo. The calculations are carried out under 
isochoric conditions, as appropriate given the above considerations. 

The thermal calculations are shown in Fig. [13] for the case of L46 = 1, and demonstrate the limitations of the assumption of 
photoionization and thermal equilibrium. The gas heats up very quickly near the quasar where the ionizing flux is high and less 
so on larger scales. Typically, the gas within the inner 50kpc heats up to 2 x 10 K over a period of ^ 10"* years (cf. equations 
A3,A4). Further heating occurs on somewhat longer timescales which can be comparable to the quasar lifetime. Unless one 
happens to intercept young quasars then the inner 50 kpc region of the halo would be heated to < 10^ K over the quasar lifetime. 
We conclude that, for moderate luminosity quasars, the gas in the central regions of the haloes can be efficiently heated to 
< 10^ K. Nevertheless, further heating is delayed and may take longer than the quasar duty-cycle or other relevant dynamical 
timescales of the problem. Assuming Xq is similar for all quasars then higher temperatures can be reached in the inner regions of 
the haloes of bright quasars since the ionization timescale is shorter and the flux higher. 
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Fig. 13. — Time dependent thermal structure of the halo after a quasar ignites at t = (and a cloud at a distance r from it is then exposed at time t + r/c) 
assuming L46(?) = &{t) where 0{f) is a step function. Also plotted are estimates for typical quasar lifetimes. Time dependent photoionization effects are 
important and prevent the gas from reaching very high temperatures (close to the Compton temperature) over the quasar lifetime. Right: Ionization fractions for 
several abundant elements as a function of time. Colors denote different ions (see legend) while different lines denote gas at different distances from the quasars: 
30kpc (solid), and 300 kpc (dashed) lines. The rapid ionization of all low level ions is clearly seen and takes place over < lO'* years. Such gas will therefore 
disappear from the line-of-sight and would become transparent. We note the general agreement between the ionization timescales obtained and the analytic 
approximation (equation A3). 



We have also calculated the time-dependent ionization structure of the halo and trace several of the ionization levels in Fig. [T3] 
As shown, all the low ionization stages usually probed by absorption line studies quickly attain negligible abundances and are 
not likely to be observed unless the quasar brightened over short timescales. For example, Mg II and H I would be immediately 
ionized and the gas is quickly heated. Time-dependent ionization deviations from a steady-state solution can be important but 
mainly for the high ionization levels whose absorption lines would appear in the X-ray region of the quasar spectrum. 

The above analysis suggests that it may be possible to treat the halos of quasars as having the same structure as galaxy halos 
since little structural changes are likely to occur over the quasar lifetime by radiation effects alone. We note that this is may not 
be the case if quasar activity is a recurring phenomenon or if they have additional means of affecting their environment (e.g., via 
jets, winds, etc.). To conclude, for the purpose of this work, steady-state photoionization calculations can be adequately used to 
describe the ionization and thermal state of the halo gas at all times. 
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